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ABSTRACT 



^ ! We have embedded the texture of the neutrino mass matrix with three families into the 



MaVaNs scenario. We take the power-law potential of the acceleron field and a typical texture of 
active neutrinos, which is derived by the D4 symmetry and predicts the maximal mixing of the 



atmospheric neutrino and the vanishing Ues- The effect of couplings among the dark fermion and 

Oh" 

. active neutrinos are studied by putting the current cosmological data and the terrestrial neutrino 



experimental data. It is found that the neutrino flavor mixings evolve as well as the neutrino 



^ ■ masses. Especially, Ue3 develops into the non-vanishing one and 6'atm deviates from the maximal 
H ; 

mixing due to couplings among the dark fermion and active neutrinos. 
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1 Introduction 



One of the most challenging questions in both cosmological and particle physics is the nature 
of the dark energy in the Universe. At the present epoch, the energy density of the Universe 
is dominated by a dark energy component, whose negative pressure causes the expansion of 
the Universe to accelerate. In order to clarify the origin of the dark energy, one has tried to 
understand the connection of the dark energy with particle physics. 

Recently, Fardon, Nelson and Weiner [1] proposed an idea of the mass varying neutrinos 
(MaVaNs), in which the neutrino couples to the dark energy. Gu, Wang and Zhang [2] also 
considered the coupling of the scalar, such as Quintessence, to the neutrinos. It should be also 
noted that the variable neutrino mass was considered at first in [3] , and was discussed for neutrino 
clouds [4]. 

The renewed MaVaNs scenario [1] has tried to make a connection between neutrinos and the 
dark energy. In this scenario, an unknown scalar field which is called "acceleron" is introduced, 
and then, the neutrino mass becomes a dynamical field. The acceleron field sits at the instanta- 
neous minimum of its potential, and the cosmic expansion only modulates this minimum through 
changes in the neutrino density. Therefore, the neutrino mass is given by the acceleron field and 
changes with the evolution of the Universe. The cosmological parameter w and the dark energy 
also evolve with the neutrino mass. Those evolutions depend on a model of the scalar potential 
strongly. Typical examples of the potential have been discussed by Peccei [5] . 

The MaVaNs scenario leads to interesting phenomenological results. The neutrino oscillations 
may be a probe of the dark energy [6, 7]. The baryogenesis [2, 8, 9], the cosmo MSW effect of 
neutrinos [10] and the solar neutrino [11, 12] have been studied in the context of this scenario. 
Cosmological discussions of the scenario are also presented [13, 14, 15, 16, 17]. The extension to 
the supersymmetry have been presented in ref. [18, 19]. 

Now one needs to construct the realistic MaVaNs scenario with three families, which is consis- 
tent with the terrestrial neutrino experimental data [20, 21, 22, 23, 24]. In order to get a realistic 
MaVaNs scenario, we have embedded a texture of the neutrino mass matrix with three families 
into the MaVaNs scenario and examined neutrino masses and flavor mixings [25]. We take a 
typical texture, which is derived by D4 symmetry and predicts the maximal flavor mixing and 
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the vanishing one in the lepton sector. The effect of couphngs among the dark fermion and active 
neutrinos are studied by putting the current cosmological data [26] and the terrestrial neutrino 
experimental data [20, 21, 22, 23, 24]. It is found that neutrino flavor mixings evolve as well as 
neutrino masses. Especially, C/e3 develops into the non- vanishing one and ^atm deviates from the 
maximal mixing due to the couplings among the dark fermion and active neutrinos. 

In section 2, we present the formulation of the MaVaNs scenario with three families, and in 
section 3, we study evolutions of neutrino masses and flavor mixings. The section 4 devotes to 
the summary. 

2 Dark Energy and Three Active Neutrino Masses 

In the MaVaNs scenario, one considers a dark energy sector consisting of an acceleron field, 0a 
and a dark fermion n. Only left-handed neutrinos are supposed to couple to the dark sector. 
There are two constraints on the scalar potential of the acceleron field (pa- The first one comes 
from the observation of the Universe, which is that the present dark energy density is about 
0.7pc, Pc being the critical density. Since the dark energy is assumed to be the sum of the energy 
densities of neutrinos and the scalar potential in the MaVaNs scenario, 

^dark = + ^i^a) , (1) 

the first constraint turns to 

pI + ^(0°) = O.Tpc , (2) 
where "0" represents a value at the present epoch, and 70% is taken for the dark energy in the 
Universe. 

The second one comes from the fundamental assumption in this scenario, which is that p^ja^]^ 

is stationary with respect to variations in the neutrino mass. This assumption is represented by 

dpv ^ dV{4)a{m^i)) ^ ^ 
d E m„i d E rut^i 

For our purpose it suffices to consider the neutrino mass as a function of the cosmic temperature 
T [5] . Since one can write down generally 



the stationary condition eq. (3) turns to 



where 



= - [ 



1 /-oo dyy'^^y^ + ^f 



(5) 



(6) 



ey + 1 

We can obtain the time evolution of neutrino masses from the relation of eq.(5). 

Since the stationary condition should be also satisfied at the present epoch, the second con- 
straint on the scalar potential is 



d(f)a 



(7) 



m^i{ct>a)=mli{<f°), T=To 



fl d^i d(f)a 

where Tq = 1.69 x 10~^eV. It is found that the gradient of the scalar potential should be negative 
and very small. 

One can also calculate the equation of state parameter w as follows: 

4 - h{T) 



w + 1 



1 + 



(8) 



where 



(9) 



In order to calculate the evolution of neutrino masses and we assume the power-low 
potential for the scalar potential T^(0a) as follows: 



(10) 



where A and k are fixed by the condition of eqs.(2) and (7) if the magnitude of the acceleron 
field at present , 0° is given. Then, we can calculate evolutions of neutrino masses and w. 

3 Neutrinos Masses and Flavor Mixings 

Let us discuss the neutrino mass matrix. Suppose only left-handed neutrinos couple to the dark 
sector, which consists of an acceleron field, 0o and a dark fermion n. A Lagrangian for the dark 
sector and the neutrino sector is given as 



L = VLarrf^on + A^^nn + VLaM^^PRfS + T^laMfC'^VR^ + h.c. , (11) 

i/L and vr are the left-handed and the right-handed neutrinos, respectively. Since we consider 
three families of active neutrinos, mass matrices m% M^'^ and Mf are 3 x 1, 3 x 3 and 3x3 
matrices, respectively. Then, the neutrino mass matrix is given as the 7x7 matrix 

/ mo Md\ 



M 



mjj X(j)a 

\Ml Mr J 



(12) 



in the {ul, n, ur) basis. We assume the right-handed Majorana mass scale Mr to be much higher 
than the Dirac neutrino mass scale Md, and Xcpa to be much higher than the scale of mo- Then, 
the effective neutrino mass matrix is approximately given as 

= MnM^'M^ + . (13) 

Furthermore, we have one sterile neutrino with the mass A0o and heavy three right-handed 

Majorana neutrinos. 

The first term in the right hand side of eq.(13) is the ordinary neutrino seesaw mass matrix, 
which is denoted by M^, and so it depends on the flavor model of neutrinos. On the other hand, 
the second term originates from couplings of the left-handed neutrinos to the dark sector. The 
matrix is parametrized as 



/ a \ 



niD — D 



(14) 



h 

\ c y 

where coefficients a, h and c are introduced to indicate relative couplings to three flavors. Then, 
the second term is written as 



, r^-, I ab ac 
A0„ " A0„ [ J 



(15) 



In this paper, we assume that the ordinary neutrino seesaw mass matrix M,^ dominates the 
effective neutrino mass matrix M^. Therefore, the contribution from the dark sector is the next- 
leading term. The first term depends on the model of the neutrino mass matrix. In order to find 
the effect of the acceleron on neutrino masses and mixings, we take the following typical texture 
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of neutrinos in the flavor basis of the charged lepton: 



jX Y 
Y W 
\Y U 



Y 

U 

w 



) 



(16) 



which leads to the neutrino flavor mixing matrix: 



U 



/ Ci2 Si2 \ 



(17) 



with arbitrary three neutrino masses rfji, m2 and rhs. Here, Ci2 and Si2 denote cos6'i2 and sin ^12, 
respectively. The mass matrix of eq.(16) was derived by the discrete symmetry D4 and examined 
phenomenologically [27]. This mass matrix gives Ues = and = — 1/-\/2 as seen in eq.(17). 
On the other hands, $12 is an arbitrary mixing angle. Therefore, this neutrino mass matrix is 
attractive one because it guarantees the one maximal flavor mixing and one vanishing flavor 
mixing. 

In this paper, we discuss the effect of couplings of the left-handed neutrinos to the dark 
fermion on the neutrino mass matrix with the flavor symmetry. Since the contribution of the 
dark fermion depends on couplings a, b, c, which denote flavor couplings, we examine typical 
three cases ofa = 6 = c = l (the flavor blind), a = l,b = c = (the first family coupling) and 
a = b = 0,c = 1 (the third family coupling). The other choices of a, b, c do not provide drastic 
changes of resits in these three cases. 

3.1 The case ofa = 6 = c= l (the flavor bhnd) 

At first, we consider the case that couplings among the left-handed neutrinos and the dark 
fermion arc the flavor blind. Therefore, we take a = 6 = c= lin eq.(14), which leads to the 
effective neutrino mass matrix 



Since the structure of the neutrino mass matrix is not changed by the dark sector as seen in 
eq.(18), the flavor mixings do not deviate from Ues — and — —1/V2. Assuming that the 
contribution of the dark sector in eq.(18) is small compared with the active neutrino mass matrix. 



1/ 



fx Y 
Y W 
\Y U 




(18) 
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we get the mass eigenvalues mi(i = 1 ~ 3) and eigenvectors Ui in the first order perturbative 
expansion as follows: 



Ui = Ui + u. 



(1) 



(19) 



where 



'(1) 



u 



(1) 



u: 



1 1 1 
1 1 1 
1 1 1 



u, 
a, 



A0a 

CiiUi + Ci2U2 + Ci3U3 , 



(20) 



with 



and 



ifii — frij 



{i ^ j) 



Cii — , 



Ui 



Cl2 \ Si2 

:i ' 1 

V2 / \ V2 



U3 = 



/ \ 

]_ 

V 73 / 



Since we get 



A0a 



(Ci2 - V2S12)' 
(Ci2 - V2Si2){V2ci2 + S12) 





(Ci2 - V2Si2){V2Ci2 + S12) \ 
(Sl2 + V2Ci2f 

0/ 



mass eigenvalues and flavor mixings are written in the first order approximation: 

mi = mi + — — (ci2 - V2si2)^ , 

]J2 

m2 = m2 + T— (Sl2 + V^Ci2)% 
A0a 

ms = ma , 



Ue2 
Ue3 



, (Cl2- V^Si2)(V2ci2 + Si2) 
■S12 + C12 z Z 



A0a 



7712 — mi 



1 

71' 



0, 



(21) 



(22) 



(23) 



(24) 



7 



where f/gs and are free from the contribution of the dark sector. 

Since the parameters rhi{i = 1 ~ 3) and S12 arc arbitrary, these are taken to be consistent 
with the experimental data in the 90% CL hmit [20, 21, 22, 24]: 

0.33 < tan^ ^^^n < 0.49 , 7.7 x 10"^ < Am,\^ < 8.8 x 10"^ eV^ , 

0.92 < sin^ 2e^tm , 1-5 x 10"^ < Am^tm < 3-4 x 10~^ eV^ . (25) 

The bound obtained by the reactor neutrinos [23] ^reactor < 12° is also taken in our study. 

Assuming the normal hierarchy of neutrino masses, we take the following typical neutrino 

masses at present epoch: 

mg = 0.05 eV , m2 = 0.01 eV , mi = 0.0045 eV . (26) 

These values lead to Am^tm = 2.4 x 10~^ eV and Am^^^ = 8.0 x 10~^ eV, which are consistent 
with the experimental values in eq.(25) ^. On the other hand, since we expect that Ue2 does not 
deviate from su so much, we take a typical value S12 = 0.5, which gives the relevant value of Ue2- 
By using these values, we examine the effect of the dark sector on neutrino masses and flavor 
mixings. Evolutions of masses and mixings depend on the parameter of the dark sector D^/A0°, 
which is unknown one. We show the numerical results for three cases of 

n2 

— = 10^^ eV , 5 X 10-^ eV , 10"^ eV , (27) 

because the effect of the dark sector is tiny in the case of D^/Xcl)^ < 10~^eV. On the other hand, 
the speed of sound becomes imaginary in the case of D^/A0° > 10~^eV, which will be discussed 
in the section 4. 

For each case of eq.(27), we solve equations in eq.(24) at present epoch, and then we can get 
mi, 7712, "rhs, Ue2, Ufjtz and C/gs- The numerical results are presented in Table 1. It is noticed that 
and C/e3 are free from the effect of the dark sector, while m2 and Ue2 have the seizable 
contribution from the dark sector, and mi has the tiny one. Therefore, Am^^j^^, Amg^^ and Ue2 

are time dependent in the universe. 

^There are other choices for m2 and mi. For example, the set of m2 = 0.009 eV and mi = 0.001 eV also gives 
Am^un = 8.0 X 10~^ eV. However, the numerical results in Table 1 do not so change. 
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On the other hand, the sterile neutrino mass is predicted for each case of eq.(27) as 1.0D\eV , 
2.0D^keV and lOD^keV, respectively, which depends on the unknown parameter D in the eV 
unit. 

Since we know 0a dependence of neutrino masses, we get the parameter A and k in the scalar 
potential V{(f)a) of eq.(lO). The conditions of eqs. (2) and (7) turn to be 



A 

k 



0.7pc - P° = 2.79 X 10-^^eV , 

m.i(</.J=mU,(,^0), T=TO ^^'/'a 



= 3n 



[(Cl2 - \/2Si2)' + (Si2 + \/2Ci2)' 



(28) 



where n° = 8.82 x lO^^^eV'^ is taken for the neutrino number density at present epoch . The 
values of k are also summarized in Table 1. 



Couplings 




777.3 (eV) 7772 (eV) 777i(eV) k 


Ue2 Ufja Ue3 


a — b — c — 1 


5 X 10'^ 
10-4 


0.05 0.0070 0.0045 9.47 x 10"^ 
0.05 0.0085 0.0045 4.73 x 10"^ 
0.05 0.0097 0.0045 9.47 x 10"^ 


0.59 -1/V2 
0.53 -1/V2 
0.51 -1/V2 


a = 1 
6 = c = 


10-3 
5 X 10-4 
10-4 


0.05 0.0098 0.0037 3.16 x 10-^ 
0.05 0.0099 0.0041 1.58 x 10"^ 
0.05 0.0010 0.0044 3.16 x 10"^ 


0.56 -1/^2 
0.53 -1/^2 
0.51 -1/V2 


a = b = 
c= 1 


10-3 
5 X 10-4 
10-4 


0.0495 0.0096 0.0044 3.16 x 10"^ 
0.0498 0.0098 0.0044 1.58 x 10"^ 
0.0499 0.0099 0.0045 3.16 x 10"^ 


0.46 -0.698 0.0006 
0.48 -0.703 0.0003 
0.49 -0.706 0.0001 



Table 1: The summary of output parameters and predictions of mixings. The input parameters 
are taken as 7713 = 0.05 eV, r?72 = 0.01 eV, r?7i = 0.0045 eV and S12 = 0.5 
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3.2 The case ofa = l,6 = c = (the first family coupling) 

Let us consider the case that flavor dependent couphngs of the left-handed neutrinos to the dark 
fermion. The case of a = 1, 6 = c = in eq.(14) leads to the effective neutrino mass matrix 



n2 / C12S12 

'(1) ^ I ^ ^2 

A0„ 



Mi} = — I C12S12 sl^ I , (29) 




which gives mass eigenvalues and flavor mixings in the flrst order approximation: 



mi = mi + irrCu , 



,2 



1712 = rn2 + T— S12 , 

7713 = rhs , (30) 

'Je2 — Si2 + 



Ma 7112 - mi 

1 

C/e3 = 0. 



Taking same input parameters in the case oi a — b — c — 1, we get numerical results as shown 
in Table 1. In this case, ms, and Ues are free from the effect of the dark sector as well as the 
case of a = b = c = 1. On the other hand, mi, m2 and Ue2 have the seizable contribution from 
the dark sector. Furthermore, we get 

A = 2.79 X 10-"eV , (31) 

which is the same one in eq.(28), and k is summarized in Table 1. 

3.3 The case ofa = 6 = 0,c=l (the third family coupling) 

We discuss another case of couphngs. Since the case of a = 6 = 0, c = 1 in eq.(14) leads to the 
effective neutrino mass matrix 



'12 ~Ci2<Si2 — S12 

2 

12512 C^2 "^IS 

-S12 C12 1 



= I -C12512 C?2 Ci2 I , (32) 
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we get mass eigenvalues and flavor mixings in the first order approximation: 

1 

mi = mi + T^T— Si2 , 

1 2 

m2 = rn2 + -—Ci2 , 

2 A0a 

1 

= ^^+2A0:' ^^^^ 
^'2«^2 



2 A0a 7712 - mi 
1 1 C?2 S?2 



"\/2 2-\/2 A0a \"^3 - "7,2 ma-mi^ 
1 D'^ f C12S12 C12S12 



2 X(j)a Krhs - rh2 rh^-rhi. 

Taking same input parameters in the case oi a — b — c — 1, we get numerical results as shown 
in Tabic 1. In this case, m^, and Ue3 also catch the effect of the dark sector as well as mi, 
m2 and t/e2- It is remarked that the non- vanishing t/es is obtained as well as the deviation from 
the maximal mixing of The parameter A and k are same one in the subsection 3.2. 

4 Evolutions of Neutrino Masses and w 

If 0° is fixed, in other words, D^/X is given, we can calculate the evolution of neutrino masses 
and w by using eqs.(5) and (8). In order to see these evolutions, we take D'^/X — leV^ in eq.(27) 
in this section. 

In the cases ofa = b = c= l and a = 1, b = c = 0, does not evolve because the 
dark sector does not contribute to the third family, while it does evolve a little in the case of 
a — b — c — 1. Therefore, we show the evolution of neutrino mass m2 versus the redshift 
z — T/Tq — 1 for three cases in Figure 1. 

As seen in Figure 1, the remarkable evolution of m2 is found in the region of z = ~ 1 for 
the case of a = 6 = c = 1. Evolutions of cases ofa = l, 6 = c= and a = 6 = 0, c=l are 
small compared with the one in the case ofa = 6 = c= l. Evolutions of w versus z are shown 
in Figure 2. These behaviors of w are similar in three cases. The value of w turns to positive 
around z = 10, which is somewhat different from the result by Peccei [5], in which the power- law 
potential was also taken in the one family model and then, w becomes positive near z = 20. 
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a=l,b = c = 
'■a = b = Q,c=l 



a=b=c= 1 
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Figure 1: Plot of m2 versus the redshift z for three cases. 




Figure 2: Plot of the equation of state parameter w versus z for three cases. 
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It was remarked [15] that the speed of sound, Cs becomes imaginary in the non-relativistic 
hmit at the present and then the Universe cease to accelerate. However, our prediction of is 
positive and around 0.1 at the present epoch if the condition D^/(A0(j) < 5 x 10~^ is satisfied. 
It means that one needs careful study of the next-leading term in the non-relativistic limit of 
neutrinos. 

In order to get the next-leading term in the non-relativistic limit, the function in eq.(6) 
is expanded as 



dyy"^ ^ 1 



ev + l ' ey + l 



+ ... 



Then, the speed of sound is given as ^ 

2 



E 

i=l 



( dm,, 



dz 



i=l 



5,. 



1^' 



where 



dz 



A dn„ 

i=l 



(34) 



(35) 



(36) 



dz 



'2 /•' 



zy + 1 



f 



^ rco dyy'^ 
-L JO eV+l 



6.47 . 



(37) 



2 fOO dyy^ 
JO ey+l 

In eq.(36), the time derivative is replaced with z derivative. The first term in the right hand 
side of eq.(36) is the leading term in the non-relativistic limit and is negative definite because of 
dmi,i/dz < and dn^^dz > 0. Since the numerator of the second term, which is the next-leading 
term, is reduced to 



4To E 



^ E.2 



5To E 



1 dm. 



(38) 



fi^iJ ' & dz) — V^=t 67 V^=t ei dz J ' 

the second term is positive definite. The magnitude of dm^i/dz is model dependent. In our 

power-low potential for V{(j)a), the magnitude of the first term becomes very small. Then, the 

second term is non- negligible one in the case of < (!?(10~^)eV. If we take the condition 

-DV(A0a) < 5 X 10-"^, the leading term is very small and then, becomes positive due to the 

second term. Actually, we checked numerically relative magnitude of each term in the right hand 

side in eq.(36). Thus, the positive is understandable in our model. 

paper including the detail derivation of and its physical implication will appear soon. 
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5 Summary 

We have embedded the texture of the neutrino mass matrix with three famihes into the MaVaNs 
scenario. We take a typical texture of active neutrinos, which is derived by the symmetry 
and predicts the maximal mixing of the atmospheric neutrino and the vanishing Ue3- The effect 
of couplings among the dark sector and active neutrinos are discussed by putting the current 
cosmological data and the terrestrial neutrino experimental data. 

It is found that flavor mixings evolve as well as neutrino masses. Especially, 1/^3 develops 
into the non- vanishing one and ^atm deviates from the maximal mixing due to the dark sector in 
the case of couphng a — b — 0, c — 1. It is also remarked that the speed of sound, could be 
positive in our model. 

Thus, the three family neutrino texture can be embedded into the MaVaNs scenario, where 
the power-law potential of the acceleron field, (pa is taken. It is noticed that our results do not 
change so much even if we take the exponential potential of 4>a- The related phenomena of our 
flavor mixing varying scenario will be discussed elsewhere. 
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